Slow light in dielectric composite materials of metal nanoparticles 
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We propose a method for slowing down light pulses by using composites doped with metal nanopar- 
ticles. The underlying mechanism is related to the saturable absorption near the plasmon resonance in a 
pump-probe regime, leading to strong dispersion of the probe refractive index and significantly reduced 
group velocities. By using the non-collinear scheme, it is possible to realize the total fractional delay of 
43 or larger values. This scheme promises simple and compact slow-light on-chip devices with tunable 
delay and THz bandwidth. 
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Recently, slow light [Q1 0] has attracted much attention be- 
cause of its fundamental significance and promising long- 
reaching applications including all-optical storage, switch- 
ing and data regeneration in telecommunications, enhance- 
ment of interferometers performance and of nonlinear op- 
tical processes. It has been observed using different physi- 
cal mechanisms such as electromagnetically induced trans- 
parency J3l0], coherent population oscillations (CPO) JD- 
7], stimulated Brillouin H] and Raman scattering |0], pho- 
tonic crystal waveguides ifiol [Till , fiber Bragg gratings 
double resonances [131 and others. The existing 
methods, however, are restricted to diverse limitations, 
most notably large propagation length (such as km-scale 
fibers for Brillouin scattering), small operation bandwidth 
(for example, 2 ~ 50 MHz in electromagnetically induced 
transparency), as well as long pulse durations and others 
that must be circumvented for a given application. In par- 
ticular, for practical applications, miniaturized designs of 
the delay line and its integration in optical networks are 
requested, that can be mass-produced and will work reli- 
ably under varying conditions. Recently progress into this 
direction has been achieved by the realization of on-chip 
delay based on different schemes lH, [l4l[l5ll . 

In the present paper we propose a new approach for the 
realization of a slow-light device which relies on compos- 
ite materials doped with metal nanoparticles (NPs). In such 
composites the absorption near the plasmon resonance be- 
comes saturated with increasing intensity lfl6j[l7ll . because 
the light-induced nonlinear change of the metal dielectric 
function leads to a shift of the plasmon resonance lfl7ll . The 
key feature of the considered system is the retarded nonlin- 
ear response of metals in the ps time range. Therefore in 
the pump-probe regime the permittivity will show large- 
amplitude oscillations, if pump and probe frequencies are 
sufficiently close. As a consequence, a THz-scale dip in- 
duces in a homogeneously broadened spectral shape, and 
(due to Kramers-Kronig relations) a steep change of the ef- 
fective refraction index near the probe wavelength, leading 
to a significant slowing down. 

We consider a dielectric composite doped with small 



metal NPs, which is illuminated by two waves: the strong 
quasi-continuous-wave pump E and the weak pulsed 
probe E pr (t) with a central frequency slightly off-set from 
that of the pump. The main feature in the optical response 
of these materials is the excitation of plasmons in the spec- 
tral range of the plasmon resonance which is responsible 
for the enhancement of the local electromagnetic field in 
the vicinity of the NPs 111 811 , The total enhanced field 
in the NPs can be written as E cnh (t) = E^e'™ * + 

I?™ 11 (t) e~ lulprt . Here we assume 
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much smaller than 



with a beat frequency Q = ui — ui 
lv and u pr . The transient nonlinear change of dielectric 
function of metal can be described as 
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where r is the electron-phonon coupling time in the range 



of 1-3 ps 119D and XW = XmK; w , -u , U) ) is the 
inherent third-order nonlinear susceptibility of the metal 
NPs at the pump wavelength. In this paper we consider 
picosecond pulses and therefore the electron-electron scat- 
tering process with a response time of few fs can be ne- 
glected. The corresponding transient dielectric function of 
the metal is therefore described by 



(t) = e m0 + x 



(3) 



pcnh 



+ 2Re [E e 



cnh* 



£»^(i + l0eK>'(ii']} 

(2) 

where e mQ is the linear dielectric function of metal. From 
the above equation, one can see that e m (t) (and there- 
fore the effective index) oscillates as cos Qt. As a re- 
sult, a part of the pump beam energy is transferred to the 
probe through two-beam coupling. Correspondingly, the 
total transmittance of the probe is increased. This principle 
shows a certain analogy with coherent population oscilla- 
tions (CPO) in two level systems, but the essential differ- 
ence here is that the two beam coupling arise by the non- 
linear excitation of surface plasmons and the saturable ab- 
sorption do not arise by population inversion but by the 
nonlinear shift of the plasmon resonance. The absorption 
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dip in this mechanism is created in nanocomposites with 
a homogeneous plasmonic absorption profile by coherent 
coupling between pump and probe, being substantially dif- 
ferent from the slow light mechanism based on spectral 
hole burning II20II . in which an absorption dip originates 
from the selective bleaching of a homogeneous line in the 
inhomogeneously broadened absorption spectrum without 
any coherent beam coupling. 

Separating the spectral components at cj and ujp r in the 
electric displacement D m (t) = e m (t) E enh (t), we find 
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For spherical metallic NPs with diameters smaller than 
10 nm, the field enhancement factor (the ratio of en- 
hanced field to the incident field) is given by x = 
2>Eh/ (s m + 2e/j), £h being the permittivity of host 
medium. Thus we can write the equation for the field en- 
hancement factor of the pump x Q (uj q ) 

3e ft 



x (u ) 



which we solve numerically lfl7h . By using the standard 
Maxwell-Garnett theory [21] and the above equations, we 
obtain the effective dielectric function for the probe. For 
larger or non-spherical NPs, we apply the discrete dipole 
approximation in combination with the effective medium 
approximation 1171] by using Eqs. ([3] 01). The term 1 + 
1+ ] Ut in Eq. ([3]) for e m (uj pr ) leads to a narrow spectral dip 
in the loss coefficient and the refraction index, the width of 
which is given by r _1 , arising from the coherent coupling 
of the probe with the pump. 

Let us first demonstrate the above described principal 
mechanism of slow light for the simplest case of spherical 
metallic NPs permitting the application of the Maxwell- 
Garnett formalism. Figure 1(a) shows an example of the 
creation of an absorption dip and the corresponding strong 
dispersion of the refractive index in a layer of silica glass 
doped with Ag nanospheres smaller than 10 nm for a pump 
with an intensity of 400 MW/cm 2 at the plasmon resonance 
wavelength of 414.7 nm. The dielectric function of sil- 
ver is taken from Ref. 12211 and its nonlinear susceptibility 
x (3)= -(6.3 - il.9) x 10" 16 m 2 V~ 2 [23]. The refractive 
index of silica is 1.47 and the electron-phonon coupling 
time r = 1.23 ps [19]. The filling factor is 3 x 10~ 2 . 
The transmittance for the probe strongly increases near 
O = and drops for higher O. As the plasmonic dephas- 
ing time is in the range of few femtoseconds, the width of 
the dip is more than 100 times narrower than the plasmon 
bandwidth. This narrow transparency window in the plas- 
monic absorption spectrum leads to strong dispersion and 
a large decrease of the group velocity near O = due to 
the Kramers-Kronig relation. Note that in certain parame- 
ter regions the imaginary part of e m (uj pr ) can be negative 
and correspondingly the absorption dip change over into a 
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FIG. 1: (Color online) Creation of absorption dip and strong dis- 
persion of refractive index (a) by frequency-dependent coherent 
energy transfer from pump to probe, maximum fractional delay 
(b) and transmission (c) as functions of the pump wavelength 
for different pump intensities in a silica glass layer doped with 
very small Ag nanospheres. In (b) and (c), the filling factor is 
2.5 x 10 -2 , the propagation length is 2 /im, t = 1.23 ps. 



gain peak. For spherical particles this occurs if the condi- 



tion Im [e m (u)q)] < 



xL 3) (ES nh ) 

purely real Xm > requiring stronger pump in practice due 
to the non-zero imaginary part of xin ■ F° r spherical NPs 
with the above given parameters, this condition requires a 
rather high intensity in the range of 10 GW/cm 2 . However, 
as will be shown later for gold nanorods the same effect 
occurs also for non-spherical particles. Since the plasmon 
resonance for nanorods can be shifted to the IR spectral 
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range exhibiting the much higher value of I Xm I , a gain 
peak is already created for even much smaller intensities in 
the range of 0.3 MW/cm 2 . 

Figures 1(b) and (c) show the fractional delay F = 
Tde\/T and the transmission T as functions of the 
pump wavelength for different pump intensities. Here 
T dcl (uj pr ) = L [rig 1 — c" 1 ) is the total delay, L the prop- 

v g = c/n g 
I diUi pr the group 



agation length, c the light velocity in vacuum, „ g 



pr ) 



the group velocity, n g = d[n (to 

index with the effective index n (oo pr ), and T is the pulse 
duration. The filling factor is 2.5 x 10~ 2 , the propaga- 
tion length is 2 /im, and the duration of the probe pulse is 
1.85 ps, which is around 1.5 times longer than the electron- 
phonon coupling time r. The other parameters are the 
same as in Fig 1(a). One can see that the fractional de- 
lay increases with increasing intensity while the transmis- 
sion shows an opposite behavior. As can be seen, for 
pump intensities below 300 MW/cm 2 the delay decreases 
with increasing pump wavelength while the transmission 
increases. 

For practical applications, it is worth to study slow-light 
performance at telecommunication wavelengths. The plas- 
mon resonance can be shifted to a large extent by tailor- 
ing the shape of the metallic NPs, in particular nanorods 
are suitable to shift the plasmon resonance into the long 
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FIG. 2: (Color online) Slow light in Ti02 film doped with Au 
nanorods with a diameter of 20 nm and a length of 66 nm for 
pump intensity of 6 MW/cm 2 at 1550 nm. Other parameters are 
the same as in Fig. 3. In (a) and (b), probe pulse evolution and 
optical delay are shown, respectively. In (b), blue dotted line is 
the incident probe pulse, green dash-dotted and red solid lines are 
probe pulses corresponding to propagation lengths of 4 and 5 fim. 



wavelength range. In this wavelength range, noble metals 
have extremely large inherent nonlinear susceptibilities, e. 
g. for gold —1.5 x 10~ n m 2 V~ 2 , which is about 4 or- 
ders larger than in the visible range l24ll . To take into ac- 
count the decrease of the pump intensity during the propa- 
gation due to the residual absorption and the change in the 
complex amplitude of probe, here we simultaneously solve 
the coupled propagation equations for the cw pump and the 
pulsed probe, using the effective dielectric function, calcu- 
lated by the method mentioned above. The time-domain 
probe pulse shape then is obtained by using the Fourier 
transformation. 

Figure 2 displays an example of the incident and de- 
layed pulses for a pump intensity of 6 MW/cm 2 at 1550 
nm. As a slow light material, we take a Ti0 2 film doped 
with Ag nanorods with a diameter of 20 nm and a length of 
66 nm, exhibiting a surface plasmon resonance at around 
1540 nm. Pump and probe beams have the same propa- 
gation direction and polarization. The pulse duration was 
1 .85 ps , corresponding to 1.5r. Figure 2(a) shows the evo- 
lution of probe pulse. Though the absorption for the pump 
is saturated, the output pump beam is attenuated down to 34 
kW/cm 2 (not shown), which is ~ 200 times smaller than 
the input pump intensity. The probe is sustained primarily 
by the coherent energy transfer from the pump and then it is 
attenuated by linear loss after the strong decrease of pump 
intensity. The probe energy can become even larger than 
the incident at a certain propagation length. The possible 
occurrence of a gain peak for the probe under certain con- 
ditions was discussed above; here it arise already for rather 



low intensities. This can be clearly seen from Fig. 2(b); the 
peak intensity of the probe after propagation over a length 
of 4 /xm (green dash-dotted line) is nearly 1 .2 times larger 
than that of the incident pulse. For a propagation length of 
5 /jm, the fractional delay is about 1 .79, corresponding to 
a delay-bandwidth product of 2.68, with an energy trans- 
mittance (the energy ratio of output pulse to the incident 
pulse) of about 0.76. As the pulse width is broadened dur- 
ing the propagation in the medium, the peak intensity was 
decreased slightly stronger than estimated from the trans- 
mittance and was about 0.55. 

The group velocity is 198 times smaller than the light 
velocity in vacuum. The pulse distortion was evaluated to 
be around 0.76 by using the formula 1 25 , \M, 
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where I m and I out are input and output probe pulse inten- 
sities normalized to have the same maximum, respectively, 
and At is the peak delay time. Because of the large in- 
herent nonlinear susceptibility of gold in this wavelength 
range, the large fractional delay is achieved even with rel- 
atively low intensities of few MW/cm 2 . 

Figure 3 presents the dependence of fractional delay F 
and transmittance T as functions of the pump intensity in 
the same medium as in Fig. 2 at 1550 nm. The figure shows 
that the fractional delay decreases with increasing pump in- 
tensity, while the transmittance increases. In particular, for 
higher pump intensity, the transmittance for the probe may 
exceed unity, indicating amplification by an energy transfer 
from the strong pump. This feature is a significant advan- 
tage compared with most of the known slow-light mecha- 
nism. Pulse distortion, in this case, varies from 0.8 to 0.66 
with increasing intensity (not shown). Fig. 3 shows that for 
a given propagation length of 5 /im a fractional delay in the 
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FIG. 3: (Color online) Dependencies of fractional delay F (blue 
solid line), transmittance T (green dotted line) on the pump in- 
tensity in TiC>2 film doped with Au nanorods at 1550 nm. The 
diameter and length of nanorods are 20 nm and 66 nm, respec- 
tively. The propagation length is 5 fim, the filling factor 5 x 10~ 2 , 
and the probe pulse duration 1.85 ps. 
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range of 2 (corresponding to a delay-bandwidth product of 
3) can be obtained with a transmittance higher than 0.5. 

In our knowledge, relatively few experiments in differ- 
ent schemes have measured relative pulses delays larger 
than these predictions. Comparing results obtained with 
chip-compatible design, we refer to Ref. [Q] with a relative 
delay of 0.8 or to Ref. ltl5ll with a fractional delay of 1.33. 
On the other hand, considerably larger fractional delays 
have been realized by using cascaded microring resonators 
using photonic wire waveguides (but with relatively small 
probe transmission) [ 14] or by using double resonances of 
Cs atoms lfl3ll . However the above given numerical exam- 
ples with a fractional delay in the range of 2 are not the 
physical limits of the method proposed in the Letter but 
determined by the chosen design. Next we will show that 
significantly larger fractional delays can be realized using 
a modified arrangement. 

The main limitation of the delay line in collinear config- 
uration considered up to now arises due to the limited prop- 
agation length caused by the attenuation of the pump. This 
can be circumvented using a transversely pumped waveg- 
uide as shown in Fig. 4(a) where the probe pulse is guided 
by the compact waveguide structure. The waveguide dis- 
persion could be neglected in comparison with the disper- 
sion arising from the NPs. The thickness of the film is 
taken to be 1 /zm, resulting in an attenuation of the pump 
intensity less than 10 %, and, therefore, the pump intensity 
can be approximated to be constant. The other parame- 
ters are the same as in Figs. 2 and 3 except the pump in- 
tensity of 0.28 MW/cm 2 . This corresponds to a pump en- 
ergy of 0.2 nJ for a pump pulse duration of 200 ps, which 
can be considered as quasi-continuous wave in compari- 
son with few ps duration of probe, and a transverse width 
of the waveguide of 5 /im. Note that the polarization di- 
rections of pump and probe have to be parallel due to the 
selective plasmon excitation for non- spherical NPs. Fig- 
ure 4(b) presents the evolution of the probe pulse normal- 
ized to the peak intensity of the incident probe pulse. In 
this figure, the peak intensity of the probe increases by a 
factor of 3.1 after propagation over a length L of 90 /mi. 
For detailed presentation, in Fig. 4(c) we show the nor- 
malized probe pulses corresponding to propagation lengths 
up to 90 fim with a equidistance of 6 fim. The figure 
shows a total fractional delay of about 43, corresponding 
to a delay-bandwidth product of 65 and a slowing down 
factor of around 270. The corresponding pulse distortion 
remains in an acceptable range with changes from 0.60 to 
0.94 with increasing propagation length. In this configura- 
tion, as the probe can be amplified, we have no any princi- 
pal limitation of the delay time, since with increasing prop- 
agation length and correspondingly larger pump energy the 
delay time can be increased significantly only limited by 
the available pump source. 

To conclude, we proposed and studied theoretically a 
slow-light mechanism based on composites doped with 
metal NPs. If two pulses with a frequency difference 
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FIG. 4: (Color online) Delay for probe with a duration of 1.85 ps 
and pump intensity of 0.28 MW/cm 2 in TiC>2 film with thickness 
of 1 /im in non-collinear configuration: (a)- configuration, (b) 
and (c) - evolution of probe pulse. Other parameters are the same 
as in Fig. 3. In (c), blue dashed line is the incident pulse and solid 
lines are the delayed pulses corresponding to propagation lengths 
L up to 90 (im with a equidistance of 6 fim from the left to the 
right side in the order. 

smaller than the electron-phonon coupling rate propagate 
through such a medium, plasmon-induced oscillations of 
the nonlinear permittivity arise, creating a strong disper- 
sion of the effective index for the picosecond probe pulses, 
and correspondingly a significantly small group velocity. 
This scheme can be applied over a broad spectral range 
from visible to infrared by tailoring the sizes and shapes 
of the metal NPs. We have shown that employing these 
composites in a collinear arrangement reduced group ve- 
locities of the probe with a slow-down factor in the range of 
200, a fractional delay up to more than 2.5 for probe pulses 
with a few ps duration can be realized as at telecommunica- 
tion wavelengths using gold nanorods as well in the optical 
range. Avoiding pump depletion by using a non-collinear 
configuration, the relative delay can be significantly in- 
creased in a regime of probe amplification. We have shown 
that a total fractional delay of 43 (delay-bandwidth product 
of 65) at telecommunication wavelength can be realized 
in a transversely pumped waveguide geometry. Since the 
probe is not attenuated but can be amplified, this value of 
the fractional delay can be further increased with increas- 
ing length of the transversely pumped waveguide. The ad- 
vantage of this configuration includes a high compactness, 
cost-effectiveness, broad range of applicable spectral range 
with tailoring the sizes and shapes of metal NPs, wide slow 
light bandwidth up to nearly THz, and large fractional de- 
lay. This could open the perspectives for the development 
of simple, compact, cheap and chip-compatible slow-light 
devices with large delay for applications in optical telecom- 
munication and other fields. 
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